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Abstract 



Oh! 

^T) [ We present a canonical formalism facilitating investigations of the dynamical Casimir 

effect by means of a response theory approach. We consider a massless scalar field confined 
inside of an arbitaray domain G{t), which undergoes small displacements for a certain 
00 '. period of time. Under rather general conditions a formula for the number of created 

particles per mode is derived. The pertubative approach reveals the occurance of two 
generic processes contributing to the particle production: the squeezing of the vacuum 



, by changing the shape and an acceleration effect due to motion af the boundaries. The 

I method is applied to the configuration of moving mirror (s). Some properties as well as 

■ the relation to local Green function methods are discussed. 

g : PACS-numbers: 12.20; 42.50; 03.70. +k; 42.65.Vh 

Q^! Keywords: 
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1 Introduction 

After the discovery of the static Casimir effect |1| (see, e.g. Refs. [0, and |^ for a review), 
the possibility of creating particles out of the vacuum by moving one of the mirrors (see, 
e.g. |Tl|]-[^) or both plates (see, e.g. |T0|) has been analyzed. Other authors calculated the 
radiation from one single mirror (see 0-0) and the backreaction on perfectly or nonperfectly 
conducting boundaries (see, e.g. ||^-|2^). Many intresting studies also have been devoted to 
the analysis of quantum vacuum radiation induced by moving dielectrics ||31]]. We discuss a 
general Hamiltonian formalism for an arbitary domain G{t) with Dirichlet boundary conditions 
that experiences small changes during a time interval (0,T): 

G{t < 0) = G{t >T) = Go. (1) 

In our derivations we imply that the time-dependent disturbances AG(t) = G{t) Go of the 
boundary can be considered as small with respect to some parameter e, i.e. AG = 0{e). 
Having introduced a proper definition of particles together with a vacuum state, we calculate 
the number of produced particles within the framework of response theory. The result will be 
applied to the special case of moving mirror(s) and the relation to results obtained by means 
of the adiabatic approach or local Green function methods will be indicated. 



2 Canonical formulation 



2.1 Equations of motion 

We consider a non-interacting real massless scalar (Klein-Gordon) field in Minkowski-space- 
time with the Dirichlet-boundary-conditions: $ = at dG{t) and □ $ = in G{t). In the 
following we shall not quantize the degrees of freedom assigned to the motion of the boundary. 
Focussing on the scalar field sector only, we consider the Lagrangian {h = c = 1 throughout) 

C = ^d^^d^^ (2) 
Expanding $ in terms of eigenfunctions fa{r,t) that fullfil /q, = at the boundary dG(t), i.e., 
mt) = ^qa{t)Ur,t) (3) 

a 

and calculating the Lagrangian by making use of the properties Eqs. (^)-(^) given in the 
appendix, we arrive at 

L= J rf^/: = ^g'-^^^'(t)g'+ga^a/3(t)g/3 + ^ga-Ma7W'^/37W?/3 (4) 

G{t) 

together with the eigenvalue equation 

Ur,t) = -nlit) Ur,t) (5) 
and the coupling matrix 

Mo.,{t)= J dV^^^Ur,t) . (6) 

Git) 

In view of the orthonormality of the fa (see Eq. (^) in the appendix; dG = dV; dG = v dA ) 
and the required boundary-conditions the Aia/s turn out to be antisymmetric: 



Ma^it) + Mpait) = I dG-[Ur,t)Uf,t)] 

G{t) 

= f^Sia,/3)- I dGUf,t)fpif,t) = (7) 

dGit) 

In Eq. (^) and in most of the following formulae we drop the summation- and integration signs 
and declare that one has to sum over all multi- indices like etc. that do not occur on both 
sides of the equation. Introducing the canonical conjugate momenta 



Pa = WT- =qa + qpMpaif) 



(8) 



the Hamiltonian takes on the form: 



Hit) = ^pl + ^nl{t)ql+p^Ma(,it)qp . (9) 

There are two effects which could lead to an unstable vacuum: The nonstationary eigenfrequen- 
cies fla{t) due to a dynamical change of the shape of the domain G{t) - we shall refere to this 
effect as "squeezing" of the vacuum - and the additional g^-p^-coupling A4ai3, indicating the 
motion of the boundaries - the "acceleration" -effect. The total energy of the $-field is given 
as the integral over the time-dependent domain G{t) of the energy density Too 

E{t)= I dVToo = \pl + \nl{t)ql (10) 

G(t) 

of the minimal coupled energy-momentum tensor: 

T^u = d,<!>d,<!>-^g,,dp<^d'<^ . (11) 

Comparison with Eq. (^ reveals that 

H{t) = E{t) + W{t) (12) 

holds. The time-dependent transformation $(r, t) — > qa{t) results in the difference between the 
Hamiltonian, describing the time-evolution of the g^, and that of $ (which is equal to the field 
energy) : 

^//[$,n,t] =E[$,n,t] =E[g„,p^,t] . (13) 
2.2 Quantization 

Now we perform the usual canonical quantization, assuming the following set 
of equal- 1 ime- commut at ion-relat ions : 

[qa{t)Ap{t)] = [P^{t),Pp{t)] = Q , (14) 

[g,(t),P;3(t)] = z5(a,/3) . (15) 

Note that, because of (^ and (0) together with the decomposition 

n(f, t) = <i>(f, t) = p^{t) Uf, t) (16) 

of the conjugate momenta these commutation relations are consistent with those between the 
fields: 

[l>(f,t),l>(r-",t)] = [fl(f,t),fl(r-",t)] = (17) 

and 

[l>(f,t),n(r',t)] =i(5(r-f') , (18) 

which are valid inside the domain G{t). Equation (|TB|) and therefore Eq. ( [T^ ) are not pointwise 
equalities (think of 9G), they have to be read as identities of L2(G')-distributions. 



3 Vacuum-definition 



3.1 Interaction- representation 

For performing the pertubation-theory we shall adopt the interaction representation. Accord- 
ingly, the time-evolution of the operators will be governed by the undisturbed energy operator 
Eq = E{t < 0) = E{t > T) defined via 

= E(t) + iy(t) = io + AE(t) + iy(t) = Eo + i^'i(t) = ^o + ^iW • (19) 

Time-dependent operators A(t) obey the equation of motion 
dA OA 

^ = ^1+3^ , (20) 
while the dynamics of any given quantum state is described by 

= -^H,(tm , (21) 

where the interaction Hamiltonian Hi is specified in Eqs. (p5D-(^) below. 

3.2 The number operator 

A proper set of particle creation and annihilation operators should be introduced with respect 
to the unperturbed Hamiltonian Hq = H{t < 0) = H{t > T) = Eq according to 

a„(t) = {2nl)-'/\nlUt) + ^Pa{t)) , (22) 

together with the static frequencies = Vtait < ^) = ^ait > T). 
They obey the equation of motion: 

— = ^£^0, OqJ = -W^aCta (23) 

and the commutation relation 

[a„(t),a+(t)] = 5(a,/3) . (24) 

These particle creation and annihilation operators diagonalize the unperturbed Hamiltionian, 
i.e., expressed in terms of the corresponding number operator 

Na = a+(t)a,(t) (25) 

it thus takes the form 

Eo = (iV« + I) . (26) 

Evidently we define the vacuum |0) as the ground state of E^: 

Va : a„ |0) = . (27) 

So Na counts the physical relevant particles in the mode a before and after the motion of the 
boundaries (due to |0,m) = \0,out)), a particle definition during the movement is not so easy 
to obtain (see section 5 and 6). 



4 Particle creation 



4.1 Response theory 



Now we investigate the change of the state vector which satisfies the initial condition 

m < 0)) = |0) , (28) 

due to a small time-bounded but otherwise arbitrary motion of the boundary by computing 
the number of created particles per mode a to first non-vanishing order perturbation theory. 
Equation (|2lD can be formally integrated by means of the time ordering operator T 



r 



exp 



dtHi{t) 



E 

n=0 



—I 







|0) 



(29) 



dtn-- - \ dtiT H^{tn) ■ ■ ■ Hl{tl) |0) 
n\ Jo Jo ^ J 



Assuming small perturbations we shall keep only the lowest-order terms of the expansion above. 
The time-evolved vacuum state reads 



l-ij dtHi{t) 





\Q) + 0{Hl) . 



In view of the property 
Va |0) = 



(30) 



(31) 



the number operator has no linear response; the first non-vanishing order is quadratic. But 
other operators A with A |0) 7^ such as e.g. the components of the energy-momentum tensor 
T^y, possess a linear response: 



(^(T)|A|^(T)) = (0|A|0)+*y dH<d\[H^{t),Am + 0{Hi) . 



To be complete we note here the general expression for the quadratic response: 



(32) 



{i^{T)\A\i^{T)) = (0|A|0)+*y rft(0|[i/i(t),A]|0) 



T T 

dt J dt' (^(0| Hi{t)AHr{t') |0) - - (0| {T[Hrit)H,{t% A} |0)) + 0{Hf) . (33) 


We are now in the position to calculate the number of particles Na created in the specific mode 
a after the time duration T, when the boundaries are again at rest. We have to evaluate the 
matrix element 



T T 



{^{T)\N^ \^{T)) = dt dt' (0| Hi{t)NaHi{t') |0) + 0{Hf) = iV„ + 0{Hl 



(34) 



with the interaction Hamihonian 
H,{t) = AE{t) + W{t) , 

where 



(35) 



(36) 



and 



W{t) = q^{t)M^p{t)p^{t) . (37) 
According to Eq. (|^) M.aa(t) = holds and thus the cross terms vanish, i.e. 

Va {0\W{t)NaAE{t')\0) = . (38) 
As a consequence, one obtains 

T T 

Na = J dtj dt' ((0| AE{t)Nc.AE{t') |0) + (0| W{t)NaW{t') |0)) . (39) 



Up to quadratic order the sqeezing- (first term) and the acceleration-effect (second term) de- 
couple, so that: = + N^. With the aid of Eqs. (|6D and (|3) we get: 



T T 



.,1 



dt j dt' -Anlit) An'^it') (0|g2(t)iV„g2(^/)|o) 


T T 

+ J dtj dt' M,x{t)M,r{t) {0\Ut)Px{t)N^qAt')Pr{t')\0) 





(40) 



Evaluation of the expectation values by utilizing the equation of motion (in the i?o-dynamic) 
leads to 



T T 



dt J dt' ^±^Anl{t) Anlit') expmllt' - t]) 
+ J dtj dt' S^p{t)S^p{t') exp(z[fi° + nl][t' - t]) 




T T 



(41) 





with the symmetric matrix: 



So^pit) = Sp^{t) = -Maf^it) 















\ 





(42) 



A more compact form of Eq. (^Tj) is derived by using the Fourier-transformation JF 
with ip(t) <pii-^) = [^^]{'-^) '■ 



= ^^,mU2Ql)\' + \S^,iQl + = + . (43) 

This spectral representation above provides the main result of the perturbative approach: The 
number of particles created in the mode a (with an energy Q^) decomposes into a squeezing and 
an acceleration contribution. The occurance of these two distinct contributions reflects the two 
basic degrees of freedom that characterize dynamical changes of the boundaries: deformations 
in shape and motion of the boundaries. The total number of particles produced after the time 
interval (0,T) is obtained by summing/integrating over all modes a: 



N = ^N^ . (44) 

a 

However, this quantity is in general ill-defined and requires an appropriate regularization. This 
may be most easily achieved by introducing explicitly a frequency cut-off which simulates im- 
perfect conducting boundary conditions. Whether a cut-off-independent contribution remains 
after its removal will depend on the particular configuration of boundaries under consideration. 
We should also note, that for "well behaving" time dependencies the spectral form (^) may 
provide sufficient convergence of the summation/integration in the limit of perfect conductors. 



4.2 Discussion 



Formula ( ^31) gives the produced number of particles in a given mode a in lowest order 
pertubation-theory in Hi. If we want to analyze a particular situation within the perturbative 
approach we need to specify the very general quantity Hi. Furthermore, we have to introduce 
the magnitude of the displacement of the boundaries e. By inspection, we see that all quanti- 
ties entering Hi are at least of order e, i.e., Af2^ = 0{e) and = 0{e) (and A/q, = 0{e) 
as well) reveals that Hi{t) = 0{e). This already indicates that both, the squeezing- and the 
acceleration-effect are of the same order of magnitude. As a matter of fact, the adiabatic 
approximation applied frequently in studies cannot be sufficient in general. In the adiabatic 
approach the acceleration term (second term) is not obtained. Consequently, the dynamics 
reduces to a set of decoupled ordinary differential equations of the form: 

x^it) + nl{t)x^it) =0 . (45) 

They can be solved by means of a scattering-theory approach (see e.g. [|l5l). The adiabatic 
approach thus mainly accounts for the squeezing contribution. We like to stress that there 
are indeed situations, where is much smaller than N^. Let fla(t) have a functional form 
corresponding to a reflectionless potential, for example: 

one obtains = 0. Or, consider e.g. harmonic oszillations of the boundaries with the fre- 
quency ujq over a long period T with Tuq ^1. In that case only particles with = ^f^o 
are produced by the sqeezing effect, but the acceleration effect also leads to the production of 



particles with other frequencies. 

This example already demonstrates the high-resonant character of the squeezing term, while 
the accereration term does not have this property due to the summation/integration over the 
modes. As mentioned above, another feature of the Fourier transformation leads to the re- 
sult, that all expressions like total number of particles, total energy, etc. are convergent after 
summation/integration over all modes, if (and only if) the time-dependend function of the 
displacement of the boundaries is smooth enough. 



5 Squeezing effects 

To investigate the pure squeezing effect it is convenient to neglect any explicit time-dependence 
of the system. In the following considerations we restrict to situations, where the time- 
dependence enters only implicitly via some global length parameter x describing the changes 
of the shape of the boundary dG. In general x could abbreviate a set of suitable parameters 
characterizing the dynamics of the shape. While the pure acceleration effect can appear sep- 
arately e.g. as a rigid motion of the boundary as a whole, there seem to be no possible G{t) 
that will lead in a pure sqeezing effect only, since there is no shape-changing without inducing 
motions of the boundaries. Accordingly, the VL^ become functions of x ^"^^ therefore also the 
creation- and annihilation-operators as well: 

a«(x) = {2n^{x)r"\^a{x)qc. + iVa) ■ (47) 

The same holds for the total energy of the field: 

E{x) = ^a{x){NM + l)=EN + Ez , (48) 

where Ez denotes the zero-point energy and that of the contained particles. 
Equivalently, 

Va : a,(x)|0(x))=0 (49) 

for the vacuum as the ground state of E{x)- Now we investigate the potential arising from the 
constrained vacuum of the quantized scalar field, which may be expanded around a fixed but 
arbitrary (shape) configuration X = Xo + ^X'- 

V{x) = (O(xo)l ^(X) |0(xo)) = Vixo) + (1^) ^X + l Ax^ + O(Ax^) . (50) 

Formally we can introduce a force via 

which should lead to the well-known Casimir force after appropriate regularizations have been 
performed However, the second- derivative term 



with 



{0(xo)l(^~) |0(xo)) = "a (O(xo)l Wxo)} (53) 



and 



c}a„ 1 dflr 



at (54) 



ax dx " 

gives rise to an additional parabolic potential due to the fact that 

(OWI0|O(X))7^^(O(X)I^(X)|O(X)) • (55) 

As a consequence, after the summation/integration over all modes a is performed the additional 
parabolic potential has an infinite strength (easy to verify e.g. for the moving mirror example; 
in contrast to the parabohc potential the Casimir force turns out to be convergent because its 
divergent parts cancel, if both sides of the mirror are taken into account) which counteracts to 
any displacements of the boundary. This leads to the conclusion that for an empty closed system 
of perfectly conducting boundaries at zero temperature it would be impossible to observe the 
static Casimir effect via measurents of the Casimir forces exerted on the boundary. As we can 
see in (^3p, any finite change of x would lead to an infinite amount of produced particles so that 
the backreaction would compensate the Casimir force even after an infinitesimal displacement. 



6 Velocity effects 

The pure acceleration effect can be studied by chosing G{t) = fj{t) + Go- (That means a time 
dependent translation of the "rigid" domain Gq; another possibility could be a rotation.) In 
view of 

W{t) |0) = '-S^p{t)ata+p |0) ^ (56) 

the ground state of Eq is not stable under the time-evolution of if, even for constant velocities. 
Therefore, for nonvanishing r/, the diagonalization of Eq does not yield a proper definition of 
creation/annihilation operators at/cLa that describe physical (i.e. Lorentz-invariant) particles. 
Only in the frame where the boundaries are (globally) at rest a reasonable definition of particles 
by diagonalization of Eq is possible. One should notice that: 

H{t) = EQ + W{t)= J c/y Q[n2 + (vl>)'] - n(^v)$) (57) 

G{t) 

and therefore: 

lH= I f/v(f^,A;(r)))^^= J dVfo.A''M= I dT.^aV , (58) 

Git\ G(t) G(t) 



where A = ^^p{ri) denotes the Lorentz transformation with the four-velocity = 7(1,7/) = 
(1, 0, 0, rf) / ^/l^^rf of the boundary. 

Also for 7 ^ 1 the Hamiltonian as defined above does not coincide with the energy operator : 
H^Eo' = jdV'{f')^^ = j dV'[K+fk)^^ = jdV'u^f^^u" (59) 

Go Go Go 

introduced by a co-moving observer. This is an indication of the fact, that H{t) describes the 
dynamics for the observer time t, and not for the time t' of the co-moving frame. 



7 Moving- mirrors- configuration 

Now we are going to apply the formalism derived above to the special case of two parallel 
mirrors placed aX z = ri{t) and z = rjlt) + l{t) or one single mirror located at z = ri{t). 
(We assume 7]{t < 0) = T]{t > T) = and l{t < 0) = l{t > T) = /q.) 
With I{t) = {r]{t);l(t) +r]{t)) it follows (see appendix): 

Aniit) = ^-^ = = {^Di m (60) 

And then the squeezing term has the simple form: 

N! = ^,m^^)\' . (61) 

For the acceleration term one does not obtain such a simple expression, an explicit form follows 
from (|S3), (!§, (|S|) and (|2p, inserted in or (||). With G{t) = I{t) = (r/(t);oo) the 
squeezing term vanishes (the same as in the case / =const) and the acceleration term yields 
after some simplifications [see (^) and (|9lD ]: 



00 

AT, = - / dnl nl ml + i^J ) f . (62) 



(|6^) enables us to calculate also the total energy radiated by a mirror 
into the l-|-l-dimensional G{t): 

E = n^N^ + 0{e') = ^Jdt m + 0{r^') . (63) 
7.1 Mechanical properties 

The force exerted upon each mirror can be calculated by computing the divergence of the 
symmetrized energy-momentum tensor: 

U = d'^f., = \id,^)m\ . (64) 



□ $ = is valid only in G{t) but it is not at dG(t). E.g., a mirror placed ai z = r]{t) with 
T]{t = 0) = and ri{t = 0) = induces the following source term (this can be verified by means 
of Fourier analysis) : 

□ $ = (nV$)5(0) , (65) 
where n denotes the normal to the plane mirror. Therefore: 

fu = d^f,, = {{d,^)nV^6{z) . (66) 

Expressions (p^)-(pB|) describe the force acting on the mirror at only one side, for a complete 
examination (and for renormalization) it is necessary to take both sides into account. 
Taking the vacuum expectation value, we obtain the mechanical force density 

f = n{0\{nV^f\0)6{z) . (67) 

The corresponding force is obtained after integration over space. So a mirror at t = experi- 
ences only the static Casimir force, also for non-vanishing i] and V. Other forces (e.g. ~?7, see 
and [^-[^) do not occur at t = but possibly at later times, when the state vector 
describing the system has changed: \ip) 7^ |0). 



8 Remaining questions 

Some modifications of the formalism presented so far become necessary, if we turn to the 
electromagnetic field (polarizations, gauge, etc.) or if Neumann boundary conditions would be 
required. However, the general structure of the formalism and the results presented remain very 
much the same. Another possible generalization of the results of this paper is to apply them to 
dynamical situations leading to different vacua, i.e. G(t < 0) 7^ G{t > T) and |0,m) 7^ \Q,out). 
As we can see in section 5 and 6, then it is nessecary to distinguish between the particle 
production due to the dynamical Casimir effect and the one which results already from the 
comparision of different vacua. The investigation of non-perfect conducting boundaries does 
not seem to be feasable in a straightforward manner within the canonical formalism presented 
above. This requires more involved studies. 



8.1 Comparision with other results 

Applying our approach to the dynamical parallel-plate configuration we also recover most of the 
results obtained earlier (see [|l^-[^]) provided the used approximations are taken into account 
carefully. E.g., for the example investigated in : G{t) = (0, L^ll+e sm{2ujit)]) for < t < T 
with ui = tc/Lq and uiT ^ 1 we reproduce the obtained result for euiT <C 1 : 

N, = N^ = ^ieuj^Ty . (68) 

For the radiation of a single mirror (see lU-ll) it is possible to compare the total radiated 
energy from (^) with the formula (3.15) of Ref. [Q: 

00 

E = {6tt)-^ [ dxo (69) 



with V = 1] and dxo = dt (The in the original formula is probably a printing error.) The 
result is the same because the additional factor of 2 in (^) is due to the fact, that (|63[) in 
contrast to ( |69D describes only the particles radiated into G{t), i.e. the emission to the right. 
It could be interesting to compare the canonical formalism presented in this paper with that of 
Ford, Vilenkin (see 0), Moore (see |11|]) and that of Fulling, Davies (see 0, and also 0). 
The main differences are: 

-In our approach the dynamics is governed by the Hamiltonian, but otherwise calculated using 
a coordinate transformation or Green's functions. 

-The regularization by means of the point splitting method, usually applied to local quantities 
as for deriving the renormalized energy-momentum tensor, is expected to provide results that 
are independent of the infinitesimal displacement vector e'^, is avoided in the canonical formal- 
ism in favour of regularization procedures (if necessary) applied on mode sums. 
Under which conditions both approaches will lead to a unique result requires careful investi- 
gations for each configuration of boundaries under consideration. E.g., one needs to clarify in 
which coordinate system should be calculated and renormalized (see section 6). 

9 Appendix 

We have to introduce a complete set of real and orthonormal eigenf unctions of the Laplace 
operator, satisfying the Dirichlet boundary conditions /q, = at dG(t): 

dVUfp = Sia,P) , (70) 



G(t) 



dV{VU){Vfp) = {n^{t)f5{a,(3) , (71) 

fain fair ') = Sir -r') m G(t) . (72) 

For the moving-mirrors-configuration the domain Git) can be expressed as follows: 

Git) = lit) ® G^ (73) 

with a time-dependent one- dimensional /(t) determining the separation and the time- 
independent subdomain G^ defining the areas of the plates. 
Then the eigenfunctions fa factorize: 

/4r,t) = /J(z,t)/,^(f^) (74) 

with a = in,r) , (3 = (m, s) and z denotes the parallel component of the position vector r. 
Accordingly, f^ = at dG-^ and /| = at dlit) . 
The frequencies decompose into two corresponding parts: 



^lit) = M' + PM 



(75) 



In G^ : 



I dV^f^f^ = 5{T,s) , (76) 
/ dV^ (V/,^) (Vf^) = {n^r S{r, s) , (77) 

/,-^(a^)/,^(6-^)=5(a^-6^) . (78) 



Note that, if we expand $ in eigenfunctions only in G-^ : 

W.t) = f^r{z,t)f^{r^) , (79) 

r 

where \/r (pr — at dl{t), the n-dimensional Lagrangian becomes a sum of effective 1- 
dimensional Lagrangians: 

L = fl dz\{4>l - {V4>rf - {n^ni) . (80) 

r I{t) 

In I{t) we may expand the 0^ according to: 

(t>r{z,t)^^q(n,r){t) fl{z,t) , (81) 
n 

that fulfil the relations 

dzfji{z,t)fl{z,t) = S{n,m) , (82) 



m 



m 

and 



/ dz{Vfl{zm^fi{z,t))^{nl{t)f5M , (83) 



^ /J(a,i)/J(6,i) = 5(a-6) in 7(i) . (84) 

n 

For instance, wc may specify the parallel-plate configuration via I{t) — {r){t); l{t)+r){t)) together 
with the eigenmodes refering to the constrained dimension: 



/«(^.*) = ^|^«m(|^l^-^Wl) (85) 
with eigenfrequencies: 

. (86) 



The coupling matrix then reads 

Ma/9W = -M/9.W = |||^a/3+||X/3 , (87) 

where ior n ^ m : 

Z,TY1T) 

X,= [(-ir+"-l]^^^^5(r,.) , (88) 



^a/3 = (-l)"^+"^^<^(r,«) (89) 



and with Qa/s = -^ap = for n = m. 

Considering one single mirror as the limiting case G{t) = {f]{t)] oo) it follows: 



where V denotes the principal value. Special care is required in any calculation involving such 
distribution-like functions (e.g. order of integration); for instance to obtain the following result: 

Ma-y{t)Mf,,{t)=fi'{t)nl5{a,P) . (91) 
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